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Abstract 

Viscous core correction models are used in free wake simulations to 
remove the infinite velocities at the vortex centreline. It will be shown 
that the assumption that these corrections converge to the Biot-Savart 
law in the far field is not correct for points near the tangent line of a 
vortex segment. Furthermore, the self-induced velocity of a vortex ring 
with a viscous core is shown to converge to the wrong value. The source of 
these errors in the model is identified and an improved model is presented 
that rectifies the errors. It results in correct values for the self-induced 
velocity of a viscous vortex ring and induced velocities that converge to 
the values predicted by the Biot-Savart law for all points in the far field. 



1 Introduction 

In engineering rotorcraft simulations such as Refs. the wake generated 

by the rotor blades is often represented by a series of connected straight vortex 
hne segments. The total velocity induced by the wake at points in the flow field is 
then found by summing the contributions of each individual segment, for which 
an analytical solution of the Biot-Savart law is readily available (such as Ref. 
More sophisticated approaches, such as Johnson (Ref. Q) represent part of the 
wake with quadrilateral panels or, as in Ref. 0] with parabolic arcs. For a point 
close to a particular line segment, the induced velocity value calculated with the 
Biot-Savart law is unrealistically high. In literature, these large velocities are 
corrected by regularising the Biot-Savart law with a model that represents the 
viscous effects inside the vortex core. One of the first researchers to introduce 
such a model was Scully (Ref. [ll[). For points located on the center line of a 
line segment, the induced velocity contribution of that segment is undefined in 
the potential case and zero for the regularised model. For a curved filament 
discretised using straight segments, this means that the self-induced velocity 
due to local curvature is not taken into account. Several methods are used in 
literature to account for this small but nonzero contribution. Scully (Ref. [ll|) 



*Submittcd for review to the Journal of the American Helicopter Society on April 17, 2010. 
Accepted for publication in the Journal of the American Helicopter Society on July 28, 2011. 



1 



uses analytical expressions for the self-induced velocity of a curved filament, 
while in other codes a cut-off method is used (Ref. [3|). A third option is 
to take no special measures for the self-induced velocity of a curved filament 
(Ref. i). 

In this note, the validity of viscous corrections as found in literature will 
be analyzed for the case of a vortex ring discretized with a set of straight line 
segments and for the case where it is represented exactly using a parametric 
curve. 



2 Filament Induced Velocity Calculations 

For a straight vortex line between point A and B, the induced velocity at a 
point rp can be calculated with (Ref. Q, P- 39), 

p /•/32 p 

V = — — / sin^d/J e = — — (cos^i - cos/32) e (1) 

where the angles /3i and /32 are called the view angles, e is the direction vector 
of the induced velocity and h is the perpendicular distance between the vortex 
line and point P (see Fig. ITa)). 

When a vortex filament can be represented with a parametric curve C(u), 
the Biot-Savart law can be written as 

^ , , F dl(u) X r(u) 

^ ' An |r(u)|3 ^ ' 

where dl{u) = dC{u)/du, the derivative of the curve at the parametric point u. 
Then, for m quadrature points on the filament the induced velocity increment 
Avi is calculated, converting the Biot-Savart integral into a finite sum. 



E 



WiAv^ (3) 



with Wi the weights of the quadrature rule. 

The efficiency and accuracy of these methods can be verified by calculating 
the induced velocity in the far field of a vortex ring and comparing the results 
with analytical solutions. For a set of points in the plane of the vortex ring on 
the X-axis (0 to 10 times the vortex radius R, see Fig.[2a|), the induced velocity 
is calculated and the relative errortQ w.r.t. the analytical results as given by 
Castles and De Leeuw (Ref. Q) for both methods are displayed in Figs. [2B1 
and[2cl 

Results calculated with the quadrature method (Fig. [2bl) for m = 32 (for a 
total of 128 Biot-Savart function evaluations per evaluation point) are accurate 
to machine precision everywhere in the field except for a small region close to 
the vortex ring itself (3/4i? 5/4i?). With this number of Biot-Savart law 
evaluations {A9 w 3°) for the same evaluation points, the segment line method 
(Fig. Ec)) only reaches an accuracy of approximately 10~^. The second-order 
trend for this method as reported in Ref. ^ is clearly visible: a ten-fold increase 
in the number of segments reduces the relative error by a factor 100. 



^The relative error is defined as |reference-approximation|/|rcfcrcnce|. 
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As the distance between evaluation points and the vortex ring is reduced 
to zero, the induced velocity increases unboundedly. On the vortex ring itself, 
the induced velocity is undefined. To remove this singularity, core models are 
introduced that mimic the viscous behaviour of a real vortex. In general, vortex 
cores have three nonzero velocity components but it is assumed that only the 
swirl velocity is non-zero. The other two components (axial and radial velocity) 
are neglected in most engineering wake models. 

On page 601 of Ref. [8|, Leishman gives a modification of Eq. ^ that incor- 
porates the Vatistas (Ref. [13]) core model, 

■p h 

- (cos /3i — cos /32 ) e (4) 



47r (r2" + /i2")i/> 

where Tc is the core size and n is a small integer. Similar modifications to the 
Biot-Savart law can be derived for other core models, Johnson (Ref. P) gives 
the modification for the model of Scully (Ref. [ll|) on page 543. Equation ([2]) 
can be modified in a similar way to remove the singular behaviour at the core 
centre, 

. . _ r dlju) X r{u) 

An (r2" -)- /i2n)l/n |r(u)|3 ^ ) 

In Fig. m the shaded regions denote the area where induced velocity values 
are modified by the two-dimensional core correction model. Near the center- 
line of an element, the modification region extends to infinity due to its two- 
dimensional nature. The actual width of the correction region depends on the 
specific core model used, for the Rankine core model it is exactly the same as 
the core size re- For other models, the width will be slightly larger than the 
core size. 

A vortex segment with end points A and B is shown in Fig. [3] together 
with thirteen paths that approach it. Along path 1 to 8, the induced velocity 
(calculated with Eq. ([T])) increases without bounds. Along path 9 the induced 
velocity is zero everywhere except for the end point, where it is undefined. As a 
consequence, at the end of paths 10 to 13 (on the vortex centerline), the induced 
velocity is zero. 

Comparison of the correction regions in Fig. [Ijwith the paths along which the 
induced velocity increases unboundedly (Fig. [3]) shows that the unconditional 
use of the perpendicular distance h in Eqs. ^ and ([5]) to correct induced velocity 
values is wrong. Corrections are applied where no asymptotes are encountered. 
A simple solution is to check where the projection of point P on line segment AB 
(P') is located (see Fig. |4a|). If it falls between the end points of the segment, 
the perpendicular distance is used as before. If the projection lies before point A 
or after point B, the perpendicular distance h is replaced by the radial distance 
to the respective end point |ri| or |r2|. In terms of the view angles /3i and /32, 
Eq. (|H) is modified as follows, 

r d 

--i^ (,g„ + rf2n)l/n ('^OS/3,-COS/32)e, 

d^lril if cos/3i<0, (6) 
where ^ d = |r2| if cos/32 > 0, 
d = h elsewhere. 
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A comparison of the shaded regions in Figs. [Ta] and |3a| shows that using 
Eq. dH]), the viscous core corrections are restricted to a region close to the vortex 
segment and do not extend to infinity. 

When the length of the vortex segment |ro| is reduced to zero, the distances 
from the end points to the evaluation point are equal, and the influence region 
reduces to a circular area (see Fig. I4bl) . Equation ^ is modified as follows, 

. ^ r |r(^)p dliu)xriu) 

47r (r2" + |r(u)|2")i/n |r(u)|3 ' ^ ^ 

The influence of these modifications on the prediction of the self-induced velocity 
of a viscous vortex ring will be assessed here. For a viscous vortex ring with a 
Rankine core model, Saffman (Ref. [lOj], p. 195) gives an analytical solution for 
the self-induced velocity, 

where Tc/R is the ratio of the core size to the vortex ring radius {tc/R -C 1). 
For a range of relative core sizes of interest to rotorcraft applications, the 
self-induced velocity is calculated using the quadrature-based method and the 
straight-line vortex method as outlined before. The segmentation results were 
calculated using a circle discretised with 360,000 segments (A0 — 0.001°) and 
the results for the quadrature-based method are generated with Gauss quadra- 
ture rules with 128 abscissae. Results of these calculations are shown in Fig. [S] 
For the original core correction model (Fig. [5a|) . it is clear that the self-induced 
velocity is underpredicted by approximately 44%, both for the segmentation 
method and the quadrature method. If Eqs. (O and Q are used (Fig. I5bl) . 
self-induced velocity values are the same as those calculated with Eq. ([8]). 



3 Conclusions 

Using a viscous vortex ring as a basis, it is shown that core correction models 
used to regularise the Biot-Savart law contain an error caused by the uncondi- 
tional use of the perpendicular distance in those formulas. Self-induced veloc- 
ities of a viscous vortex ring are shown to converge to the wrong value. The 
deviation is approximately 44% from the theoretical value, both for the classi- 
cal approach where the vortex ring is discretised using straight segments and a 
numerical quadrature method that uses an exact representation of the vortex 
ring. A simple modification to the regularised Biot-Savart formula is proposed 
for both methods that rectifies these errors. 
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Evaluation points (segm.): 
[0 ^ 3/4/?] U [5/4i? ^ lOiJ] 



Evaluation points (quad.): ■^^^ 
^ = ^ [0 ^ lOfi] n 

(a) Evaluation points on the X-axis where 
the induced velocity is calculated. 
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(b) logj^Q of relative error of induced veloc- 
ity for the parametric vortex ring method 
for Gauss-Legendre quadrature rules with 
m abscissae. 
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(c) L2-iiorm of relative error of induced 
velocity for segmented vortex ring method 
as a function of the discretisation. 



Figure 2: Convergence of induced velocity calculations of quadrature-based 



method |(b)| and segmentation method (c) to potential reference values for eval- 
uation points in the plane of a vortex ring (a) on the positive X-axis, from (0, 0) 
to (10i?,0). Analytical values (Ref. are used as reference. 




Figure 3: Approaching a vortex segment AB along different paths results in 
asymptotic behaviour for paths 1 to 8. Due to symmetry, along path 9, the 
induced velocity is zero everywhere except for the end point and as a results, 
along paths 10 to 13 no asymptotes are encountered when approaching the 
centerline. 



7 




■7- 

S 1 
if 

c 0.4 





Analytical reference values 

Quadrature-based results 

Segmentation results x 





























0.001 0.111 
Relative eore size, fv/i? 



(a) Self-induced velocity prediction with 
original core correction model (Eqs. 
and (O). 
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(b) Self-induced velocity prediction with 
modified core correction model (Eqs. ^ 
and (0). 



Figure 5: Self-induced velocity of a viscous vortex ring as a function of relative 



core size for the original (a) and modified |(b)| core correction model. 
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